Quantum Fluctuation-Induced Uniaxial and Biaxial Spin Nematics 
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It is shown that zero point quantum fluctuations (ZPQFs) completely lift the accidental continuous 
degeneracy that is found in mean field analysis of quantum spin nematic phases of hyperfine spin 
2 cold atoms. The result is two distinct ground states which have higher symmetries: a uniaxial 
spin nematic and a biaxial spin nematic with dihedral symmetry Dihi. There is a novel first order 
quantum phase transition between the two phases as atomic scattering lengths are varied. We find 
that the ground state of atoms should be a uniaxial spin nematic. We note that the energy 
barrier between the phases could be observable in dynamical experiments. 



Since the successful observation of magnetic ordering 
in condensates of sodium atoms in optical traps [ijjthe 
subject of spinor gases has attracted much attention [3,0, 
0, H S 01 • The impressive efforts to measure spin depen- 
dent interaction energies or two-body scattering lengths 
of cold atoms [1, [^, [3, 11, [l2| and recent progress in 
cooling hyperfine spin F = 3 chromium atoms [13| have 
further stimulated extensive interest in spin correlated 
many-body states of cold atoms. 

In addition, in recent work, a much better under- 
standing of spin collective phenomena of cold atoms in 
optical lattices has been achieved. Particularly among 
many other states in optical lattices[l3|, various quan- 
tum spin nematics have received much attention. Spin 
nematics have been suggested as ground states of F = 1 
sodium atoms (^^Na), F — 2 rubidium atoms 
and F = 3 chromium atoms (^^Cr) in optical lattices 



isl . 3, [13, [3 [3] ; they are also suggested to be relevant 
to = 3/2 cold atoms[2^, [2lj. Quantum spin nemat- 
ics are fascinating because of the rich topology of their 
ground state manifolds which gives rise to exotic spin 
defects and unconventional vortices with non-integer cir- 
culation integrals. Spin nematics haven't been success- 
fully realized in conventional solid state systems but it is 
widely believed that at least some of them are likely to 
be created and observed in cold atomic vapors. The spin 
nematics of bosonic cold atoms that have been discussed 
so far do not break the (lattice) translational symmetry, 
but do break rotational symmetry. The absence of lattice 
translational symmetry breaking can be attributed to the 
relatively weak super exchange interaction between two 
bosons if they are spatially antisymmetric. 

In fact, our mean field calculations show that for cold 
atoms such as rubidium (*^i?&) in the F = 2 manifold, a 
one-parameter family of spin nematic ground states, in- 
cluding ones which are characterized by completely differ- 
ent symmetries, are degenerate in energy (see discussions 
after Eq.(l3|)). In particular, uniaxial nematics that are 
invariant under any rotation around an easy axis are de- 
generate with biaxial nematics which arc invariant only 
under a dihedral group ( either Dih2 or Dih^ depending 
on eigenvalues of the nematic matrix, see below). The 
origin of this degeneracy is accidental due to the spe- 



cial form of energy in the dilute limit, very analogous 
to the highly degenerate family of solutions for d-wave 
pairing obtained by Mermin in the 1970's[^. In this 
Letter we shall study the effect of zero point quantum 
fluctuations (ZPQFs) on nematic order in cold atoms. 
ZPQFs are known to lead to a wide range of physical phe- 
nomena, including the Lamb shift in atomic physics |23|. 
the Coleman- Weinberg mechanism of spontaneous sym- 
metry breaking|2^, fluctuation-induced first order tran- 
sitions in liquid crystals and superconductor s 1251 . and 
order from disorder in magnetic systems (26i \2l\ . For 
F = 2 cold atoms, we show that ZPQFs completely lift 
the accidental continuous degeneracy appearing in pre- 
vious mean field calculations and result in two distinct 
states with higher symmetries: 

a) the uniaxial spin nematic with Z2-symmetry; 

b) the biaxial spin nematic with Dih^ symmetry. 

We consider F = 2 atoms in optical lattices using the 
following Hamiltonian: 
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Here k is the lattice site index and < kl > are the nearest 
neighbor sites, fi is the chemical potential and t]^ is the 
one-particle hopping amplitude, ip'^ is a traceless sym- 
metric tensor operator that has been introduced in a pre- 
vious work 16)]; aL,bL,CL are three on-site coupling con- 
stants which are determined from three two-body scatter- 
ing lengths ap, F = 0,2,4 and on-site orbitals. Compo- 
nents ipa/sy P — x,y,z are linear superpositions of five 
spin-2 creation operators, V'mj?i '^-F = 0,±1,±2. The 
number operator p, the dimer or singlet pair creation 
operator 'D\ the total spin operator F^ are defined as 
p = l/2trV'V, 'D^ = 1/V40tr^t^t^ = -ita.M^l^i,^^ . 

The amplitude of a condensate is a traceless symmetric 



tensor, Xa/3 =< i/'. 



>. In this condensate, atoms oc- 



cupy a one-particle spin state \x > that is defined as a lin- 
ear superposition of five F = 2 hyperfine states |2, > 
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Here Ua, a — x,y,z are components of a unit vector 
n{9,(j)); Hx — sin^^cos0, = sin 6* sin and Uz = cos 6. 
Y2,mp, rnp = 0, ±1, ±2 are five spherical harmonics with 
1 = 2. The mean field energy per site as a function of x 
can be obtained as 



E 



ah 



MF 



Tr(x*x)Tr(x*x) + -^Tr(x*X*)Tr(xx) 



Tr(x*x), 



(3) 



where z is the coordination number. Minimization of E 
with respect to the tensor x yields nematic, cyclic and fer- 
romagnetic phases studied previously. Particularly when 
cl < 0, cl < 46^ (as for rubidium atoms i^'^Rb)), min- 
imization of energy in Eq.Q requires that x be a real 
symmetric tensor (up to a phase). An arbitrary solution 
X can be obtained by applying an 5*0(3) rotation and a 
U{\) gauge transformation to a real diagonal matrix x, 



X 



/4Me*"^7^x7^~^ 



(4) 



TZ is an 5*0(3) rotation matrix, and x is a normalized 
real diagonal traceless matrix [ig, 



X = 




Tr(xx) = 2' 



(5) 



Finally, the value of chemical potential /i is 

^iL=^i + 2ztL = [aL + cl) M, M = iTr(x*x)- 

M =< pk > is the average number of atoms per lattice 
site. The shifted chemical potential fiL = fJ- + ^zt^ mea- 
sured from the bottom of the band depends on scattering 
lengths through its dependence on ol, c^. 

In the mean field approximation, all of the quantum 
spin nematics specified by different diagonal matrices 
Eq.([5]) have the same energy, i.e. are exactly degener- 
ate. For instance, when 2x2;a: = ^Xyy = —Xzz = l/\/3, 
up to an overall 50(3) rotation, all atoms are con- 
densed in the hyperfine state |2,0 >. This choice of 
X represents a uniaxial spin nematic. When Xxx ~ 
—Xyy = 1/2 and Xzz = 0, the atoms are condensed in 
the state ^ (|2, 2 > +\2, —2 >). In general, consider the 
parametrization of the matrix elements: (^ S [0, 27r] ) 



sin(C - f ) 
73 



' Xyy 



sin(g-f ) 
V3 ■ 



sin($ - f ) 
V3 



To study ZPQFs, we first examine the energy spectra 
of collective modes. For this, we expand V'^ about the 
mean field Xi 



(7) 



where the superscript ly = x,y,z,t,p labels zero point 
motions of x along five orthogonal directions: three 
50(3) spin modes (x-,y-,0-mode) for rotations about the 
X, y and z axes, respectively, a spin mode (t-mode) for the 
motion along the unit circle of and a phase mode {p- 
mode) describing fluctuations of the condensate's overall 
phase. The corresponding matrices are 




1 

1 )(8) 




and = 2x(C + f ), LP = 2x(C)- These matrices are 
multually orthogonal, Tr(L^i'') = 2(5^^. The five modes 
are decoupled and the corresponding operators obey the 



9^ 
1,1^ 



klOfj.1,. 



bosonic commutation relations. 

We expand the Hamiltonian in Eq. ([1]) using Eq. ([7]) 
and keep the lowest order nonvanishing terms. The result 
is a Hamiltonian for the fluctuations which is bilinear 
in 9j.^, 9k,u- It can be diagonalized by a Bogoliubov 
transformation. The result can be be expressed in terms 
of Boguliubov operators, 6*^ ^ and 6^ ^, 



n = 



E 



eq{2mBNvl + eq) 



(9) 
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Here eq = 4ii^^(l — cos^adi) is the kinetic en- 
ergy of an atom with crystal quasi-momentum q = 
{(lx,<ly,<lz)] dL is the lattice constant, m^jv = ^/^tLd\ 
is the effective band mass, v^, {v = x,y, z,t,p) is the 
sound velocity of the i^-mode in the small- |q| limit; 

AMtLdl {bLG^^^ - cl), vj = AMtLdl (-cl), 
AMtLd\{aL + Cl)- The velocities of three spin 
modes depend on a x-dependent 3x3 symmetric ma- 
trix G"'3, = -(l/2)Tr([L",LP] [lP.Lp]). More ex- 
plicitly, C^^ is a diagonal matrix with elements: G^^ = 
4sin2(^-^), = Asin^{i+^-), G^^ =4sin2f. No- 
tice that the velocity of the phase mode Vp can be written 
as yj HL/mBN, only depending on the band mass niBN 
and the chemical potential /i^; it is independent of ^. 

Unlike the mean field energy Em f that is independent 
of ^, the zero point energy of spin nematics per lattice 
site E is ^-dependent ( a summed over z). 



E{^) 



2Nt ^ 



^^{2vfiBNvl{i)^e^. (10) 



Here Nt is the number of lattice sites. The main con- 
tribution to the ^-dependence of energy is from fluc- 
tuations of wave vector |q| ^ ruBNVaiO'i for conden- 
(6) sates, this characteristic momentum is much smaller 
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FIG. 1: (Color online) The zero point energy of spin nematics 
per lattice site as a function of ^ (in units of \Mcl\^^'^ /t^j^'^ , 
cl < 0). Inset is g{x, y), the amplitude of ~ f ) ~ ^(0) ^ 
a function of x{= bL/cL)- y{= |cl l/fi) has been set to be 8 x 
10~^. A quantum first order phase transition occurs at &l = 
as a result of ZPQFs. For a positive 6l, the energy minima 
correspond to uniaxial spin nematics while the maxima to 
biaxial spin nematics with Di/14 symmetries. 

than h/dL. The ^-dependent energy is proportional to 
Sa"a/(^i^i)' amplitude of it can be expressed as 
IMclI^^"^ /t^/'^ g{bL/cL, Icil/iL), where gix,y) is a dimen- 
sionless function that is studied numerically. In Fig. (TJ 
we plot the energy as a function of ^. We find that there 
are two distinct phases when the parameters occur in 
the region cl < and > CL/i. They are separated by 
6l = line. When is positive, ^ = 0, 7r/3, 27r/3, ... 
are ground states. These are uniaxial nematic states 
which are rotationally invariant along an easy axis (see 
Fig. ([2] a)). When 6l is negative, the points of ^ = 
7r/6, 7r/2, 57r/6, ... are the stable ground states while the 
uniaxial nematics are unstable. 

The energy as a function of ^ has following symmetries 

E{0=E{-0,E{^+O^E{0, (11) 

These are a result of the rotation and gauge invariance 
of the energy function. They are found by examining 
the mean field solutions in Eq. (j!]) , ([6]) . A rotation TZ of 
120° around the X = y — z line effectively transforms a 
solution at ^ to ^ + 7r/3, or 7^'^x(0^= -x(^ + f )• And 
because the energy is invariant under 50(3) rotations 
and is an even function of x, one finds that E{^) = E{£_ + 
7r/3). In addition, a rotation of 180° around the x + 
y = line in the xy-plane transforms x(C) to x(~C)? 
TZ^xiOT^ = x(~C)) so similarly one finds that E{£^) — 
i?(— ^). These symmetries only depend on rotational and 
gauge invariance and are exact. So the energy is an even 
and periodical function of ^, with the period equal to 
7r/3. For an analytical function with these symmetries, 
^ = 0, 7r/6 and 7r/3 are always extrema. Our calculations 
of the zero point energy are consistent with this. The 
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(a) (b) 

FIG. 2: (Color online) (a) Wavefunctions of spin nematics at 
various ^; (b) wavefunctions of uniaxial spin nematics (U), bi- 
axial spin nematics with Dih4 symmetries (B), ferromagnetic 
condensates (F) and cyclic condensates (C). Plotted here are 
spin wavefunctions ^'s (represented by the usual spherical 
harmonics) in spherical coordinates {p,9,(j)). p = {"^i s{0, 
and *s = J2af3 ^/^Xa/snarif). ria, a = x,y,z are compo- 
nents of a unit vector n{6,(j)){see also Eq[2]); colors indicate 
the phase of wavefunctions. The phase boundaries between U 
and C, C and F, F and B are obtained by minimizing Eq.(|3]) 
(similar to the calculations in Ref.[l3]); the phase boundary 
between U and B is obtained by taking into account ZPQFs. 

above observation also indicates that up to an SO{3) 
rotation and a phase factor, states at ^ and at ^ -I- 7r/3 
are equivalent. 

We now examine the topology of the manifold of 
ZPQF-induced spin nematics. Without loss of generality, 
we first consider the uniaxial nematic state at ^ = 0, i.e., 
Xxx = Xyy = -1/2V3 and Xzz = 1/a/3; in this case, the 
nematic easy axis specified by a unit vector e is pointing 
along the z-direction (See Fig. [l](a)). Such a uniaxial 
state is invariant under an arbitrary rotation around the 
nematic axis e; it is further invariant under an inversion 
of the nematic axis e — > — e which is evident in the plot 
for the wavefunction. The vacuum manifold for the uni- 
axial nematic is therefore simply S"^ jZ^ ® S^; here S"^ jZi 
is where the nematic director lives and is the unit cir- 
cle of condensate phase variable. Unlike in uniaxial spin 
nematics of = 1 atoms [i^, here the spin orientation 
and condensate phase are not entangled. 

For the biaxial nematic state at ^ — 7r/2, Xxx = 
—Xyy = 1/2 and Xzz = 0. The state is invariant under 
the eight element dihedral-iam group Dih^. The seven 
rotations that leave the state invariant (up to a phase 
shift) are 90°, 180° and 270° rotations about the z-axis, 
180° rotations about the x- and y-axes and about the 
X ± y — lines in the xy-plane. Four of these rotations. 
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the 90°, 270° rotations around z-axis, and 180° rotation 
around x±y = lines must be accompanied by a shift of 
the phase of the condensate by tt. The manifold there- 
fore is [SO{3) X S^]/Dih4, where the Dihedral elements 
in the denominator contain the rotations mentioned plus 
the TT-phase shifts. These nematics contain half-vortices. 

On the other hand, a generic biaxial spin nematic (with 
no accidental symmetries) is only invariant under a rota- 
tion of 180° around x-, or y- or z-axis and the invariant 
subgroup is the Klein- four or dihedral- two group [iij; it 
has lower symmetries than either of the spin nematics 
selected by ZPQFs. 

The perturbative calculation carried out in this Letter 
is valid when fluctuations of the order parameter x s-re 
small. Here we estimate the fluctuations along each of 
the five directions v = x,y, z, t,p. Using Eq. [71 we evalu- 
ate the amplitude of local fluctuations in an orthogonal 
mode, Ai{i') = {0\6l ^^6 JO) . Taking into account the 
expression for v^, we have the following estimate for the 
relative amplitude of fluctuations, 

M ^ tL ' ' ^ ' 

is the effective scattering length of each mode which 
is a linear combination of scattering lengths in _F = 
0,2,4 channels: ha = ^— y — (02 — 14) — \{O'0^o,^)^ 
at = f (02 - 04) - I (ao - ^4), flp = \{2a2 + bai) + 
^ (flo — aji). For ^'^ Kb atoms in optical lattices, ha/ap ~ 
10~^, ~ 50nfc and when ^ SOOnfc, the relative am- 
plitude of fluctuations in x-, y-, z-, or t-mode is typically 
less than one percent. Finally, we would like to point 
out that in nematic insulating states the properties spin 
modes are very similar to what are discussed here. 

In conclusion, we have found that ZPQFs lift a contin- 
uous degeneracy in spin nematics. Only uniaxial spin ne- 
matics or biaxial nematics with dihedral-four symmetries 
are selected as the true ground states. For Rb^'' atoms in 
the hyperflne spin-2 manifold, b^ is between — Sc^ and 
— lOci. According to the above analysis, the ground 
state should be a uniaxial spin nematic. The ZPQF- 
induced energy landscape can be experimentally mapped 
out by investigating the macroscopic quantum dynam- 
ics of condensates prepared in certain initial states. A 
condensate of rubidium atoms initially prepared at state 
^(|2, 2 > +|2, -2 >) (corresponding to ^ = 7r/2 point), 
because of the ZPQF-induced potential shown in FiglU 
could evolve towards a condensate with atoms at state 
i|2,0 > +^(|2,2 > +|2,-2 >) (corresponding to 
^ = 7r/3 point). This leads to a temporal oscillation 
of the population of atoms at state |2, >, which can be 



studied in experiments. This study could yield informa- 
tion on the potential. We thank I. Bloch, K. Madison 
and E. Demler for discussions. This work is in part sup- 
ported by the office of the Dean of Science, University of 
British Columbia, NSERC (Canada), Canadian Institute 
for Advanced Research, and the Alfred P. Sloan founda- 
tion. 

Note added: After the submission of the manuscript, 
we notice that similar results are obtained in Ref . l30l . 



[1] J. Stenger et al, Nature 396, 345 (1998); D.M.Stamper- 
Kurn et al, Phys. Rev. Lett. 80,2027 (1998). 

[2] T. L. Ho, Phys. Rev. Lett. 81, 742 (1998). 

[3] T. Ohmi, K. Machida, J. Phys. Soc. Jpn. 67, 1822 (1998). 

[4] C. V. Ciobanu et al, Phys. Rev. A 61, 033607 (2000). 

[5] M. Koashi, M. Ueda, Phys. Rev. Lett 84, 1066 (2000). 

[6] L. Santos, T. Pfau, Phys. Rev. Lett. 96, 190404 (2006). 

[7] R. Diener, T. L. Ho, Phys. Rev. Lett. 96, 190405 (2006). 

[8] J. L. Roberts et al, Phys. Rev. Lett. 81, 5109 (1998); N. 
Klausen et al, Phys. Rev. A 64, 053602 (2001). 

[9] E. G. M. van Kempen et al, Phys. Rev. Lett. 88, 093201 
(2002). 

[10] H. Schmaljohann et al, Phys. Rev. Lett. 92, 040402 
(2004). 

[11] A. Widera et al. New J. Phys. 8, 152 (2006); A. Widera 
et al, Phys. Rev. Lett. 95, 190405 (2005). 

[12] M. S. Chang et al. Nature Phys. 1, 111 (2005). 

[13] A. Griesmaier et al, Phys. Rev. Lett. 94, 160401 (2005). 

[14] M. Greiner et al. Nature 415, 39 (2002). 

[15] E. Demler, F. Zhou, Phys. Rev. Lett. 88, 163001 (2002). 

[16] F. Zhou, G. W. Semenoff, Phys. Rev. Lett. 97, 180411 
(2006). 

[17] R. Barnett et al, Phys. Rev. Lett. 97, 180412 (2006). 
[18] J.-S. Bernier et al, cond- mat /06 1 2044 ; J. -S. Bernier et 

al, Phys. Rev. B. 74, 155 124 (2006). 
[19] K. Eckert et aZ., 'cond-mat/0603273 

[20] C. J. Wu et al, Phys. Rev. Lett. 91, 186402 (2003); C.J. 

Wu, Mod. Phys. Lett. B 20, 1707 (2006). 
[21] H. T. Tu et al, Phys. Rev. B74, 174404 (06); see also G. 

M. Zhang and L. Yu, cond-mat/0507158 
[22] D. Mermin, Phys. Rev. A9, 868 (1974). 
[23] W. E. Lamb, R. C. Retherford, Phys. Rev. 72, 241 

(1947). 

[24] S. Coleman, E. Weinberg, Phys. Rev. D 7, 1888 (1973). 
[25] B. L Halperin et al, Phys. Rev. Lett. 32, 292 (1974). 
[26] Chris. L. Henley, Phys. Rev. Lett. 62, 2056 (1989). 
[27] E. F. Shender, Sov. Phys. JETP 56, 178 (1982). 
[28] F. Zhou, Phys. Rev. Lett. 87, 080401 (2001); F. Zhou 

and M. Snoek, Ann. Phys. 308, 692 (2003). 
[29] G. W. Semenofl, Fei Zhou, Phys. Rev. Lett. 98, 

100401(2007). 

[30] A. Turner, R. Barnett, E. Demler and A. Vishwanath, 
cond-mat/0702057, 



